Anomalous Bias Dependence of Spin Torque in Magnetic Tunnel Junctions 



O 

o 

(N 
o 

Q 



loannis Theodonis,^' ^ Nicholas Kioussis/'Q Alan Kalitsov,^ Mairbek Chshiev,'^ and W. H. Butler'^ 

'Department of Physics, California State University, Northridge, CA 91330-8268 
'^Department of Physics, National Technical University, GR-15773, Zografou, Athens, Greece 
'^MINT Center, University of Alabama, P. O. Box 870209, Tuscaloosa, AL, USA 

(Dated: February 6, 2008) 

We predict an anomalous bias dependence of the spin transfer torque parallel to interface, Ty , in 
magnetic tunnel junctions (MTJ), which can be selectively tuned by the exchange splitting. It may 
exhibit a sign reversal without a corresponding sign reversal of the bias or even a quadratic bias 
dependence. We demonstrate that the underlying mechanism is the interplay of spin currents for 
the ferromagnetic (antiferromagnetic) configurations, which vary linearly (quadratically) with bias, 
respectively, due to the symmetric (asymmetric) nature of the barrier. The spin transfer torque 
perpendicular to interface exhibits a quadratic bias dependence. 

PACS numbers: 85.75.-d, 72.10.-d, 72.25.-b, 73.40.Gk 
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Theoretical calculations predict that when a spin- 
polarized current passes through a magnetic multilayer 
structure, whether spin valve[l'| or magnetic tunnel junc- 
tion, (MTj)[i, g, it can transfer spin angular momen- 
tum from one ferromagnetic electrode to another, and 
hence exert a torque on the magnetic moments of the 
electrodes. At sufficiently high current densities, this 
spin transfer can stimulate spin- wave excitations 0, Q 
and even reverse the magnetization of an individual 
domainQ. Current-induced magnetic switching (CIMS) 
has now been confirmed in numerous experiments both 
in spin valves d, 0] and more recently, in MTJs[l, Q. 
Thus, CIMS provides a powerful new tool for the study 
of spin transport in magnetic nanostrutures. In addi- 
tion, it offers the intriguing possibility of manipulating 
high-density nonvolatile magnetic-device elements, such 
as magnetoresistive random access memory (MRAM), 
without applying cumbersome magnetic fields 10]. 

While the fundamental physics underlying the spin 
transfer torque (STT) i n sp in valves has been extensively 
studied theoretically jl|, [13, [H, El, its role in MTJs re- 
mains an unexplored area thus far, except for the pi- 
oneering work of Slonczewski 0, who employed the 
free-electron model in the low bias regime. One of the 
most pressing needs is a comprehensive understanding 
of the bias dependence of the STT in MTJs, which will 
be important for the development of MRAM that uses 
CIMS for writing the magnetic memory cell. 

In this Letter, we present for the first time a compre- 
hensive study of the effect of bias on the spin torques, 
parallel [T\\) and perpendicular (Tj.) to the interface, 
in MTJs, using tight-binding (TB) calculations and 
the non-equilibrium Keldysh formalism. We predict an 
anomalous bias dependence of the spin torque, contrary 
to the general consensus. We demonstrate first that de- 
pending on the exchange splitting, T\\ may exhibit an 
unusual non-monotonic bias dependence: it may change 
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FIG. 1: Schematic structure of the MTJ, consisting of left and 
right semi-infinite FM leads separated by a thin nonmagnetic 
insulating system containing N atomic layers. The magneti- 
zation M2 of the right FM lead is along the z-axis, whereas 
the magnetization Mi of the left lead is rotated by angle 6 
around the j/-axis with respect to M2. 



sign without a sign reversal in bias or current, and it 
may even have a quadratic bias dependence. Second, 
by generalizing the equivalent circuit in Ref. Q using 
angular-dependent resistances, we show that Ty satisfies 
an expression involving the difference in spin currents 
between the ferromagnetic (FM) and antiferromagnetic 
(AF) configurations. This expression is general and in- 
dependent of the details of the electronic structure. Our 
numerical results both for the TB model and the free- 
electron model (not presented here) confirm the validity 
of this relation for any parameter set and bias. Third, 
the spin current for the FM (AF) alignment is shown to 
have a linear (quadratic) bias dependence, whose origin 
lies in the symmetric (asymmetric) nature of the barrier. 
The interplay of the spin currents for the FM and AF 
configurations is the key underlying mechanism that can 
lead to a rich behavior of the STT on bias. Finally, we 
find that the bias dependence of T± is quadratic. 

The MTJ under consideration depicted in Fig.l con- 
sists of a left and right semi-infinite noncoUinear FM 
leads, separated by a nonmagnetic insulating (I) spacer 
containing N atomic layers. The right FM lead is mag- 
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netized along the z axis (M2) of the coordinate system, 
shown in the inset of Fig. 1. The magnetization of the 
left FM lead (Mi) lies in the x — z interfacial plane, 
i.e. it is rotated by angle 9 around the axis y (nor- 
mal to the FM/I interfaces) with respect to M2. The 
chemical potentials of the right and left leads are shifted 
by the external bias, eV = fJ.h — Mr, where the charge 
current is positive when it flows along the y axis from 
left to right. The Hamiltonian for each region is de- 
scribed by a single orbital simple-cubic TB model with 
a nearest-neighbor (NN) spin-independent hopping term, 
ta, and a spin-dependent on-site energy term, e^, where 
a= L, R, and I, refer to the left, right, and insulat- 
ing regions, respectively. In the present calculations the 
left and right FM leads are identical with an exchange- 
splitting, Al(r) = £[(p) - £[(p). We use Aj = 0, 

4(R) - Ef = 1.2cV, ei- Ef = 5.4eV, TV = 5, and 
'•R(L) = ti = ^R(L),i — t = 0.4 eV, where the Fermi en- 
ergy Ep = QeV, and the tL,i and ^rj are the NN hopping 
matrix elements at the two FM/I interfaces [13] • Note, 
that Al(r-) , refers to the local effective s-d exchange inter- 
action, with values between 0.2 eV to 2 eV [l^ depending 
on the material. Under applied bias, Er ~ = the 
potential inside the insulator, ei^„ = £1 
linearly with layer number n. 

The one-electron Schrodinger equation in spin space 
for each uncoupled region a is(l6j 
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where p and q are atomic sites indices in region a, ek|| is 
the energy of the in- plane wave vector, k||, of the Bloch 

state, is the spin-dependent retarded Green's func- 
tion for each region, and / is the 2x2 unit matrix. The 
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[(el + £a) ^pq + ta {Sp.q+1 + Sp^g-i)], 

and SHpq = |Aq,, describe the spin-average and the spin- 
split part of the Hamiltonian, respectively. 

Having determined the gp^'^ for each uncoupled sub- 
system from Eq. ([T|), one can calculate the retarded 
Green's function for the entire coupled system, by solv- 
ing a system of Dyson equations which couples the gp^'^ 
through the hopping matrix elements t^ j and ti j^ at the 
two interfaces. In order to calculate the non-equilibrium 
Keldysh Greens function, G^i^f we have extended the 
approach of Caroli et al. [13| in spin space using 2x2 
Greens function matrices. 

The charge current density is [l^ 
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and the spin current density between sites i and i-f 1 is 
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FIG. 2: (Color Online) Bias dependence of the current- 
induced perpendicular component of the net spin torque per 
unit area, T±(y) — T±{0), for 6 — n/2, and various values of 
T±(0) is related to the exchange coupling energy between 
left and right FM leads. Inset: bias dependence of T±{V). 



where cr ~ (cTj., cTj,, ct^) is a vector of the Pauli matrices. 
Both and are conserved across the MTJ, while the 
spin current is not conserved (V • I^^^ ^ 0), due to spin- 
dependent scattering caused by the local exchange field 
inside the FM leads . Conservation of the total an- 
gular momentum implies that the spin current lost at an 
atomic site is transferred to its local magnetic moment, 
thereby exerting a local STT jl^, HH on site i in the 
right FM lead per unit area defined by 
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where the second equality represents the discrete form of 
the divergence of the spin current. The local spin torque, 
Ti, shown in Fig. 1, has components parallel (Ti.||) and 
perpendicular {Ti^±) to the interfacial plane. The z- 

component of vanishes because /l"' , ,^ ^ = l''/^ — 



A(/T_/i). BothT,, 



(z,i-t-l),2: (i— l,i),2: 

2^^-^ ^ /■ j^"''" and Ti^± oscillate with different 
phase and decay with distance from the I/FM interface, 
as in the case of spin valves [ill EBl ■ 

The net STT transverse to the magnetization on the 
right FM lead, a quantity which presumably experiment 
measures, is the sum of local torques, T^, 

oc 

^ ^ 2^\-^i-i,i ^ -1,0 -'00,00 -1,0' 

1=0 

where the subscripts -1 and refer to the last site in- 
side the barrier and the first site in the right FM lead, 
respectively. In the above equation, /^''oo = because 

(s) 

the components of transverse to M2 decay to zero 

as i ^ 00 [ll[. Thus, the net spin torque exerted on the 
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FIG. 3: (Color Online) Bias dependence of the parallel com- 
ponent of the net spin transfer torque per unit area, T\ | , for 
9 = 7r/2, and various values of e''^. The curves (symbols) refer 
to the left- (right-) hand ordinate, respectively. 



FIG. 4: Equivalent circuit for spin-channel currents, with 
angular-dependent resistances. 



quantized along M2, where multiple reflections within 
the barrier are neglected 18]. Substituting the currents 
^L(R) ™ I^iS- 4 into Eq. (5) of Ref. 3, we obtain 
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M2 X (Ml X M2), (6) 



right FM lead is simply the spin current calculated at the 
I/FM interface [11]. 

In Fig. [2] we show the bias-dependence of the current- 
induced perpendicular component, T±{V) — T±{V = 0), 
of the net spin torque for 9 = 7r/2, and various values of 
eK The inset displays the bias dependence of Tl(F). We 
find that T±(V) varies quadratically with bias, as orig- 
inally suggested, but not calculated, by Slonczewski Q. 
The equilibrium {V = 0) value of T±(V), related to the 
interlayer exchange coupling energy decreases in ab- 
solute value as the exchange splitting A is reduced. 

In Fig. [3] we display the bias dependence of the parallel 
component of the net spin torque, T||, (curves associated 
with the left ordinate) for = 7r/2 and for the same val- 
ues of as in Fig. [21 The most striking and surprising 
feature of T[| is its non-monotonic bias dependence, which 
can vary from almost linear to purely quadratic behavior, 
depending on the exchange splitting A. The quadratic 
bias dependence of T|| for = 1.5 eV persists even for 
small bias. Interestingly, for ^ 2 eV and 2.25 eV, T|| 
reverses its sign without a sign reversal in bias or cur- 
rent. This anomalous bias behavior may have important 
practical implications, since it suggests that the CIMS in 
MTJ may not require reversal of the current. Note, that 
T± and Tn, are comparable in size in MTJ, in contrast 
to metallic spin valves, where T± ^T^^ [ll| . 

In order to understand the underlying mechanism re- 
sponsible for the bias-dependence of Tji , we have gener- 
alized the equivalent circuit for MTJ using angular- 
dependent resistances, R^^^iO) = R'^{0)cos'^{9/2) and 
R'^'^^ie) = R''{TT)sin-^{e/2), as displayed in Fig. H The 
angular dependence of R"^'"^ (6) is equal to the inverse 
probability [P"'-'^ (^)]~^ for an electron with spin state 
\a > quantized along Mi to tunnel to a spin state \a' > 



where &\Tr) = ^{IHtt) - I^{tt)) and li'\Q) = 
^(/^(O) — are the spin-current densities for the 

AF and FM configurations, respectively. This important 
result is quite general, independent of the details of the 
electronic structure, and reduces the calculation of T^{6) 
simply to the evaluation of the spin-current densities for 
the FM and AF configurations [3]. The angular depen- 
dence of both T± and T| | is proportional to sin 6 [2] , in 
contrast to metallic spin valves [l|,[Ii|. 

Defining the dynamic current polarization [l^, P{6) — 

[iT{e) - n{e)]/i{e), where i{e) = i\e) + n{e) is the 

total current, Eq. (O reduces to 



h P{tt)I{tt) ~ P{0)I{Q) 
Ye 2 



M2 X (Ml X M2). (7) 



In order to confirm Eq. (0) or Eq. ([7]) we display also in 
Fig.[3]the bias dependence of ^[P(7r)J(7r) - P(0)/(0)]/2 
(symbols associated with the right ordinate) for the same 
values of e-'-, where the agreement is excellent. 

In order to elucidate the atomistic origin of the bias 
dependence of T\\ in Eq. ([S]), we plot in Fig. O the spin- 
current densities, I'f^ijx) and /i*''(0), versus bias for the 
AF and FM orientations, respectively for — 2.0 eV. 
One can clearly see that {Q) {I^z\i^)) varies linearly 
(quadratic) with bias for V < 0.5 eV. For the FM 
and AF configurations the tunneling can be considered 
as the superposition of two independent spin channels. 
This different bias behavior can be understood on the 
basis of the tunnel model [20| for asymmetric barri- 
ers, generalized so as to take into account both spin 
channels. The bias dependence of 7°' can be written 
as ]20], /'"(V) = /i($'")V-/2(l>'^)A$'"V2+0(V3), where 
/i and /2 are functions of the average barrier height. 
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FIG. 5: (Color Online) Bias dependence of the spin-current 
density, li''\6), for the FM and AF orientations, respectively. 
In the FM case the majority and minority electrons tunnel 
through a symmetric barrier (lower inset) with different bar- 
rier heights, , from the bottom of the band. In the AF case 
(upper inset), the two spin channels tunnel through asymmet- 
ric barriers with the same average barrier height, <1>'^, but with 
a barrier asymmetry, A^"^, of opposite sign. 



asymmetry. Here, "J*2(2) spin-dependent barrier 

height at the left (right) interface. 

In the FM configuration, the majority and minority 
electrons tunnel through a symmetric barrier (lower inset 
in Fig. [5]) but with different barrier heights, <i>'^, for each 
spin channel. In this case, ^ $i and = = Q. 
Thus, both P (0) and I^{Q) vary linearly with V, and 
hence also lif\o). On the other hand, in the AF config- 
uration (upper inset in Fig. [5]), both spin channels tun- 
nel through asymmetric barriers with the same average 

barrier height, <I>^ = 'i'^ but with barrier asymmetry of 
opposite sign, = — A<I>^. Hence, the linear bias de- 
pendence of (tt)—I^{tt) vanishes identically, and /i'''' (tt) 
exhibits a quadratic bias dependence. 

Thus, the interplay between the linear and quadratic 



bias-dependence of /z (0) and Iz'{n), respectively, 
in Eq. ^ is responsible for the non-monotonic bias- 
dependence of T|| in Fig. [3l This competition can be 
selectively tuned by varying A, giving rise to a wide range 
of rich bias behavior. For example, the purely quadratic 
bias behavior for = 1.5 eV arises from the fact that 
/^(O) — /^(O). Contrary to the free-electron model [1], 
this behavior is possible within the TB model due to the 
bell- like form of the density of states at the interfaces [3 ■ 

In summary, we predict an anomalous bias behavior 
of Til in MTJ, which varies with the exchange splitting 
in the FM leads. Tn may exhibit a sign reversal with- 
out a corresponding sign reversal of the bias, or even an 
unexpected quadratic bias dependence. The underlying 
mechanism for the unusual bias dependence is the inter- 
play between the bias dependence of the spin currents 
for the FM and AF configurations. The origin for the 
linear (quadratic) bias dependence of the spin currents is 
the symmetric (asymmetric) nature of the tunnel barrier 
for the FM (AF) orientations. We should emphasize that 
the non-monotonic bias behavior is not associated to the 
simple TB model; other systems with more complex elec- 
tronic structures can also show this behavior, provided 
that the condition /i''''(0) < Iz^\it) is satisfied. On the 
other hand, T±_ exhibits a quadratic bias dependence. 

An experimental test of our prediction can be achieved 
by relating the (observable) critical voltage for switch- 
ing measured as a function of external magnetic field to 
the spin-torque. This can be done for a magnetic ele- 
ment that has been characterized as to damping factor, 
anisotropy, and magnetic moment. Future work will be 
aimed to include the results of these calculations as an 
input into the Landau-Lifshitz-Gilbert equation, to cal- 
culate the critical current for the CIMS. 
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